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Abstract. It has been found the exact solutions for nonstationary distribution of the temperature 

■ in the liquid ring with two viscosities {p and //o) and two free boundaries of the ring. 

o' 

We consider the problem about rotational symmetric movement of viscous and nonviscous liquid in a 
ring for nonclassical model of hydrodynamics (see [HH]). A similar problem for classical models Navier- 
£> ■ Stokes and Stokes was considered in papers [3J [H 03 [51 [7] ■ Previous papers contain the proves for theorems 

of global resolvability and uniqueness, and contain the asymptotes for solutions behavior. In the moment 

■ t = 0, we see that the liquid fills a ring R 2 q < r < R w and has distribution of velocities: 

u r = ipo r^ 1 , ue = vo(f). 

Here, i?io, -R20 > 0, and tpo are constants, voir) is an arbitrary smooth function. Since ug is not depend 
from 9, we see that the specified association u r from r is unique. 

Let i?io, -R20 are f ree boundaries. In this case the liquid ring begins movement by inertia. Since we 
have rotational symmetry of an initial condition, we may assume that this one is similarly in all time 
agreement. Using this statement, we get that free boundaries are remained circles during all time of 
movement: r = Ri(t), r = Rait). 

Let us note the equations of movement and thermal conduction for nonclassical liquid [H [2] with the 
supposition of its incompressibility and homogeneity: 

+ (u • V)w MAu + -Vp = 

IT) i dt p p 

div u = 

Ln ■ dr 

pC p — = kAT+pvH 

Here, u is a vector of velocity with components p is a hydrostatic (equilibrium) pressure; p is a density; 
H is a dissipativc function; C p is a thermal capacity; k is a factor of thermal conduction; p = p v is a 
\ viscosity; M is a tensor of viscosity: 

M=( * Mo 

Po is nondissipative viscosity which can accept values with any signs. The stress tensor: 

T = —pi + 2MB 

D is the tensor with components: 

1 / dui duj 



Let r = Ri (t) be an exterior boundary of the ring and r — R2 (t) be an interior boundary of the ring. 
Since, 

u r = u(r, i); ug = u(r, t); p = p(r, t) , 
then a components of the stress tensor in polar coordinates are: 

du ( dv v 

T rr = ~P+ 2^^- + MO q 

or \or r 

(dv v\ du 

T r e = P I -j. 2^0-^-; 

\or r J or 

m , o u (dv v 

Tee = ~P+ %P Po -5 

r V or r 



Using this representation, we obtain: 
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Here u is the radial component of velocity; v is the azimuthal component of velocity. Since the stress 
vector is equal to zero on free boundaries, we get boundary conditions: 



T — 



P + 



2i/$(t) 
d (v 



■"olfr-z) = °> ifr = i2i, 2 (t); 



^ = ^(^-2^-12=0, ifr = J2 1>a (t). 
Also, we have dynamic condition on free boundaries: 

dRi(t) *(t) 



dt 



Ri(t) 



, (i = 1,2) 



Supposing /ip + /i 2 ^ 0, we have: 
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$(0) = $ ; 

u(r,0) = voir); 



The equality 



i? 2 (i) - ,R 2 (i) - i? 2 - #20; 
is the conservation law of square. The equality 



Ri(t) 
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Rio 



Rio 



r 2 v(r,t)dr= / r 2 v(r 1 0) dr 



r vo(r) dr; 
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is the conservation law of impulse moment. 

Let us write down the equation for temperature T: 



pC p 



dT dT v dT 
dt dr r d9 



= k 



1 d ( dT\ 1 d 2 T 



r dr \ dr J 
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T| t=o = T o (r,0). 

Boundary conditions for T will be specified additionally. 
Let 

■"-20 -"-20 -"-20 

v $ R\ 

u= ; v=— r ruj; * = -; a = — 5- - 1 

e = ^ ; e = T 

Jo 

be a new coordinates then we have a problem: 

duo + 2* ^ _ <9 2 w +8 3w_ 

<9r £ + 77 9?y 2 977 ' 



w l T =o = w o(»?), < jj < a; 

^ + e (eT^ = ' if?? = ' a; T - 0; 

dcj a*(* - 4) 1 } ( -> . duj\ , 

» ' w 2 + As— drj; 



dr £(f + a)ln(l + f) ln(l + f)7 V d V 
^ = 2*; *(0) = #„; C(0) = 1; £=-: 



fci?? 6 /96 „ ,d 2 Q 1 s 2 e 



4/J 3 V £>V + ^ + ??) ^r; 2 + 4 (£ + 77) 96 12 

+ (eT^ + (e + ??) V 

It is easy to see that 

u = T ^-\ cjofo) = <5 = * = 4; £(t)=8t + I 

C + 77 1 + 7; 

is an exact solution of the problem ((!])- ©. Using ©,([7]), and 6 = Q(r(t), r](r, t)), we get 

o 9 fro ^ d ®\ 16 (1 + s 2 ) 



A = ■ — 7, , D = ■ — = , £ = 



dr drj \ drj J (8T + 77+I) 

Here 

4/j 2 ' 1 ~ v 

are dimensionless parameters. 

Using methods of Group Analysis (see [9j [TOl [TT], [12]), we obtain infinitesimal operator 

3 ,j 9 1 d 
ot orj 9B 

with coefficients: 

e = C 1 t 2 /2 + C 2 t + C 3 ; 
e = C\T (4r + v + 1) + C 2 (V + 1) - 8 C 3 ; 
77 1 = (C 4 - d (r + (A/i?)77) /2) 6 + 6 2 (r, 77); 



Here Ci are arbitrary real constants (i — 1,2,3,4). A function bi{r, rj) should be a solution of following 
equation: 

16(l + e 2 ) / fry „ . Ci / A . , 

1 (C 2 + d) - — r - -r? (12) 



(8r + 77 + l) 2 V ' 2 V' 

Supposing that Cj = 0, C 2 = 0, C 4 = 0, and 6 2 (t, tj) = 0, we get: 

X = C 3 — 8 C3 — 

ar or] 

h = 8 r + rj + 1 

z 2 = e 

Using the theorem for representation of invariant solutions through invariants from |5], we receive for 
unknown function O: 

e = i 2 = <p{h) = v{8t + v + i) 

Combining this with (jHJ) we get the ordinary differential equation: 



jMji) , m o^ difih) , 16(l + £ 2 ) 



Solving this equation, and checking solution by input equation we obtain the invariant solution of equation 

9 ^>" c '- (87TfriyCT (13) 

Let Ci = 0, C?2 = 1, C4 = —2, and function 62(1", ^) is the partial solution (TT3")) : 

w ^ r> 16(l + e 2 ) 

^)=^- (8r + 7? + 1)(jB + 8A) 

then for infinitesimal operator (|10p (and for invariants), we obtain: 

X = (r + C 3 ) + (1 + t? - 8 C 3 ) ^ + (6 2 (r, r?) - 2G) A 

J! = (l + r7-8C 3 )/(r + C 3 ) 
■7 2 - (r + C 3 ) ((t + C 3 )e ~ tC 5 /2) - tC 3 C 5 /2 + 
+ 16 t(t + C 3 )(l + e 2 )/((8r + 77 + 1){8A + B)) 

Using Ji, J2, we have (in system of computer calculus Maple) the invariant solution for equation (|8j): 

fA{l + v -8C 3 )-B{ T + C 3 y 



Q{r,rj)=K 



V (r + C 3 ) 3 

A(l + 77-8C 3 )\ /8r + 77 + l 



X6XP ' B(r + C 3 ) J V r + C 3 

+ ^ 16 ( 1 + g2 ) (14) 

2 (8t + V + 1)(B + 8A)' { ' 

B (8A + B) 

Here K\, K2 are constants. It is easy to see that the expressions (113)) . (fill are solutions for equation ©, 
and: 

T — *- + GO Z 



Let 



be a boundary conditions. 

Using (fl4|) and (fT5)l we receive: 



e 1 (r) = 6(T ) a) 

e 2 (r) = e(r,o) 



0i (r) =K 
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16(l+e 2 ) 



) exp {- 
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\ / 8T+Q+1 

J V ^ 



(8 T+a+l)(B+8A) 
0„M — K ( A(l~SC a )-B (r+C 3 ) 

u 2 (t)-k y (T+ c 3 ) a 

, C5 _ 16(l+e 2 ) 
" l " 2 



) exp{- 



A(1-8C 3 ) 1 / 8r + l 



B(t+C 3 ) J \^T+C 3 



(8t + 1)(B+8A) 

The difference of temperatures on boundaries at r = is denned by the constant: 



(15) 



(16) 



C = e l( 0) - 0,(0) = jm^k + # exp {-^^j x 

(A(l + a-8C3)-BC 3 )exp{-^} (^) ¥ - (A(l - 8 C 3 ) - S Cs) (^)" 



We have distribution of temperature 0(r, 77) for r = 0: 

9(0,17) =K p'y- ) exp {- A(1+ B " C - 8C3) } + 



C5 _ 16(l+£ 2 
1 2 (tj+1)(B+8A) - 

Let us choose constants as follows: 



(17) 



C = 0, A=|, B = 6, E =i o=l, C 3 = ^ 

— C 4 5 
X = J— 3 , r = --^-3 « -0.00027127 (18) 

4(8a3-l)e X p{l-j^}-(16C3-l)exp{l-^} 6 83 

Values of temperature on boundaries at r = with such choice of constants are identical. 

Considering values of constants (|18p and expressions for temperature on boundaries (|16p . we have the 
following problem: 

90 d ( . ,dO\ 80 . , 

-5---^- 8(8r + 7? + l)— )- 1 r 7Z TT7 = ° 19 

cV di] V 977/ 3(8t + ?7+1) 2 v ; 



<-n-) 
0// 



_ 5(l + 16r) 

5 \ , ( 5 \ A-(8r + l)(8r+3) \ ^ 
3(8r+2)^ y + \^6(8r+2)a ) \ (87+Tp J 



rj—a 



C5 / 5 \ / , 2 \ / \ . ^ 



@ M=f -uj r - l)exp( -" )+ (^ri)j (21) 

Substituting values of constants (fl"8"|) in expression (|14p . we receive the solution of the problem (fTTH) - 
(ED): 

V ' " 2 V6/ V(8r + r/+l) V (8r + l) 4 / P \ 8r+ljy V ' 

Supposing temperature nonnegative, we get restriction on value of a free constant: C5 > 5/3. 
Considering 0, we get the solution (TTJI in (t, r) co-ordinates: 

T(< ' r) = Toi "l J x 



x exp 



8A 

~B B(vt + C 3 R 2 ) j \{ut + C 3 R 2 ) 

>2 
'20 



C 5 T 16T R 2 (l + e 2 ) 



Also, we have: 



2 r 2 (B + 8A) ' 

B% (t) =8vt + Rj , R 2 10 (t) = 8i4 + R 2 



10) 



4v Av 
r ' r 

As a result, exact solutions are received, the general view of Lie group for the equation (19) is received, 
the asymptote is found, and profiles for one of the solutions are constructed. These results allows to build 
another solutions for a problem about temperature distribution in the liquid ring. 
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